1 Mount Hood Environmental, PO Box 4282, McCall, Idaho, 83638, USA
2 Biomark, Inc., 705 South 8th Street, Boise, Idaho, 83702, USA
3 Washington Department of Fish and Wildife, Fish Program, Science Division, 1111 Washington Street NE, Olympia, Washington, 98501, USA
4 Mount Hood Environmental, 39085 Pioneer Boulevard #100 Mezzanine, Sandy, Oregon, 97055, USA

Correspondence: Michael W. Ackerman <>

Keywords: Northern Pikeminnow; Chinook Salmon; predation; mark-recapture; bioenergetics

Abstract

Predation of emigrating juvenile salmonids by piscivorous fishes is perhaps an under-appreciated source of mortality hindering population recovery in the Columbia River basin. Reservoirs and tailraces associated with dams on the Snake and Columbia rivers have created favorable conditions for fishes known to consume juvenile salmonids, including the Northern Pikeminnow Ptychocheilus oregonsis. Indeed, Northern Pikeminnow have been estimated to consume 8% of the approximately 200 million juvenile salmonids emigration through those rivers, annually. Less is known about the interaction between upstream, river-dwelling Northern Pikeminnow and their impacts on salmonid recovery above these reservoir systems. In this study, we examine predation by Northern Pikeminnow on juvenile salmonids, focusing on Chinook Salmon Oncorhynchus tshawytscha, in a slow-water reach of the Salmon River, Idaho known as Deadwater Slough. The reach contains historically important overwinter rearing and spring emigration habitat for Chinook Salmon and other ESA-listed species, and was formed by the failure of a small mining dam in a reservoir in the late 19th century. We estimate the abundance of Northern Pikeminnow in Deadwater Slough during peak emigrations of Chinook salmon, and further, use a bioenergetics approach to estimate the consumption potential of the Northern Pikeminnow population to assess potential impacts to Chinook Salmon emigrants. Finally, we evaulate potential impacts to adult returns of local populations. We found that… (insert concluding findings here. It is likely that predation of juvenile Chinook Salmon in Deadwater Slough has a consequential impact on ESA-listed Chinook Salmon populations in the Upper Salmon River MPG hindering local recovery efforts.

Introduction

The Snake and Columbia rivers drain a network of inland tributaries that provide essential spawning and rearing habitat for anadromous Pacific salmon Oncorhynchus spp. and steelhead O. mykiss. These fishes were historically abundant throughout the basin and consequently have tremendous ecological, cultural, and economic value. Over the last century, salmonid stocks throughout the Snake and Columbia rivers have undergone significant declines related to habitat modification in the freshwater environment (e.g., removal of beavers, river channel simplification, water withdrawals, hydrosystem development) and poor marine conditions, both affecting their survival and recruitment at multiple life stages. Consequently, action agencies have made considerable investments in the rehabilitation of tributary ecosystems and evaulation of factors attributed to the species decline (Roni et al. 2018).

One potentially important, but perhaps under-appreciated source of salmonid mortality is predation on emigrating juveniles by piscivorous fishes. In the lower mainstem Snake and Columbia rivers it is estimated that predation on out-migrating salmonids during peak emigration has a significant negative impact on the overall population and success of population recovery (Fresh et al. 2003; Winther et al. 2020). Dams and reservoirs in the Columbia River are the primary locations associated with high rates of piscine predation on salmonids (Petersen 1994; Ward et al. 1995). There are generally two mechanisms that explain these high predation zones. First, migration (movement) rates of juvenile salmon are reduced during reservoir passage (Venditti et al. 2000), thereby increasing the time migrating smolts are vulnerable to predation. Second, reservoirs and downstream tailraces associated with dams on the Snake and Columbia rivers create favorable slow-water habitat for predatory fishes known to consume juvenile salmonids, such as the Northern Pikeminnow Ptychocheilus oregonsis. Northern Pikeminnow tolerate and thrive in relatively warm, slow water habitats (Wydoski and Whitney 2003) and consequently have benefited from dams on the Columbia River, becoming abundant predators of salmonid outmigrants (Knutsen and Ward 1999). Indeed, Northern Pikeminnow are estimated to consume 8% (16.4 million) of the approximately 200 million juvenile salmonids emigrating through the lower Snake and Columbia Rivers, annually (Beamesderfer et al. 1996).

Most predation studies have focused on reservoirs (Murphy et al. 2021) and mainstem reaches (Tabor et al. 1993; Ward et al. 1995; Shively et al. 1996; Zimmerman and Ward 1999) of the Columbia River basin whereas considerably less is known about the interaction between upstream, river-dwelling piscivorous fishes and their impacts on salmonid recovery above the reservoir systems (Rubenson et al. 2020). Upstream habitats containing slower water velocities and other attributes that support piscivorous predators may overlap with essential habitat for some salmonid species and life stages. One such example is Deadwater Slough, an approximately 1.9 km long reach of unnaturally slow and deep water in the Salmon River, Idaho. The Deadwater Slough is located in a river reach containing historically important overwinter rearing and spring emigration habitat for juvenile Chinook Salmon O. tshawytscha and is inhabited by piscivorous predators, including Northern Pikeminnow and Smallmouth Bass Micropterus dolomieu. Further, the slough is part of the migratory pathway of the endangered Snake River Sockeye Salmon O. nerka population (Axel et al. 2015) and several populations of threatened Snake River steelhead.

Quantifying mortality associated with piscine predation in upstream habitats of the Salmon River is of particular interest for spring/summer-run Chinook Salmon (hereafter Chinook Salmon) in the Upper Salmon major population group (MPG). The Salmon River was historically the most productive tributary for Chinook Salmon in the Columbia River basin (Nemeth and Kiefer 1999), but those populations have become depleted in recent decades. Currently, the Upper Salmon MPG supports eight independent, extant populations, including Salmon River (above Redfish Lake Creek), Valley Creek, Yankee Fork Salmon River, East Fork Salmon River, Salmon River (mainstem below Redfish Lake Creek), Pahsimeroi River, Lemhi River, and North Fork Salmon River (NOAA 2017). The Upper Salmon MPG are stream-type Chinook Salmon and exhibit two distinct migration tactics; downstream rearing (DSR) and natal reach rearing (NRR) (Copeland et al. 2014). The DSR migrants leave their natal area as subyearlings between June and November and typically overwinter in downstream, mainstem habitats until the following spring when they emigrate to the ocean as smolts. Alternatively, NRR migrants remain in their natal areas for approximately one year after emergence until emigration to the ocean as smolts. Diversity of migration tactics provides a mechanism for coping with adverse conditions in freshwater rearing and migration environments and buffers against catastrophic events, thereby increasing population resiliency. Because the Deadwater Slough supports rearing for DSR migrants and is part of the migration corridor for NRR migrants, and juveniles from all Upper Salmon populations must rear in or migrate through it, it represents a significant habitat for the Upper Salmon MPG.

Recent studies have examined the downstream movement, distribution, and (apparent) survival of juvenile salmonid emigrants through the Salmon River, including the Deadwater Slough reach. Sockeye Salmon migrating through the mainstem Salmon River during spring were estimated to have substantially reduced survival in the Deadwater Slough reach compared to adjacent reaches (Axel et al. 2015). Similarly, emigrating DSR Chinook Salmon had a 10% reduction in transition probability through the Deadwater Slough compared to surrounding reaches during the fall and early winter months (Ackerman et al. 2018; Porter et al. 2019). Low survival in these cases was attributed to low-velocity water delaying movement rates, thereby increasing predation risk.

In this study, we estimated the abundance of Northern Pikeminnow in the Deadwater Slough and evaluated potential impacts to juvenile salmon emigrants, focusing on DSR and NRR Chinook Salmon. We hypothesize that high densities of piscivorous predators in the Deadwater Slough may explain the reduced survival (or apparent survival) observed for juvenile Chinook Salmon and Sockeye Salmon. To test this, our objectives were three-fold:

  1. Estimate the abundance (or relative abundance) of potential predators in the Deadwater Slough during the peaks of fall (DSR) and spring (NRR) juvenile emigrations.
  2. Document predation on juvenile Chinook Salmon during the DSR and NRR emigration periods using gastric lavage.
  3. Use a bioenergetics approach to estimate consumption potential (grams) of the Northern Pikeminnow population at Deadwater Slough to assess their potential impacts to DSR and NRR Chinook Salmon emigrants.

We also consider how consumption of juvenile Chinook Salmon emigrants at Deadwater Slough by Northern Pikeminnow may impact adult returns to the Upper Salmon MPG.

Methods

Study Site

The Deadwater Slough is an approximately 1.9 km section of the mainstem Salmon River located approximately 6 river km downstream from the town of North Fork, Idaho (Figure 1). The downstream end of the slough is located at the confluence of Dump Creek and the Salmon River. Around 1897, the failure of a small mining diversion reservoir in the Dump Creek drainage resulted in an erosion event that deposited substantial amounts of sediment at the confluence of the Salmon River (Emerson 1973) forming an alluvial fan. The fan created a slow, deep section in the river, spanning approximately 30 acres (12 hectares), averaging 68 m width, and resembling a small, narrow pond or reservoir.

Northern Pikeminnow Population Size

We estimated the population size of piscine predators in the Deadwater Slough using a mark-recapture survey and a catch per unit effort (CPUE) approach. Predators were sampled near the peaks of the juvenile DSR and NRR emigrations. Our initial intent was to sample during the emigration of the 2018 brood year (fall 2019 and spring 2020). However, due to logistical concerns during the onset of the Covid-19 pandemic, the spring 2020 effort was delayed until fall 2020. Although this provided the benefit of an additional abundance estimate during the DSR emigration, estimates during the NRR emigration were still lacking, and so an additional effort was added in spring 2021. Fall efforts included sampling across a two-week span, but due to budget constraints the spring 2021 effort was limited to a single week. Fall efforts were constrained within a two-week window to minimize predator emigration/immigration and to constrain the population to meet assumptions of the mark-recapture models.

Multiple capture methods were employed during the initial fall 2019 effort to reduce selectivity and bias for species and size classes. Methods included raft electrofishing, fyke netting, snorkeling, and angling. Ineffective methods were then abandoned in subsequent years. After evaluating all methods, angling proved to be the most effective method for capturing piscine predators while also minimizing potential impacts to ESA-listed adult steelhead present in the Salmon River during our sampling efforts. Therefore, the following analyses will focus on fish captured by angling, unless otherwise noted. During the fall 2019 effort it additionally became clear that Northern Pikeminnow were by far the most prevalent piscine predator in Deadwater Slough and so is our focus.

Our study relied heavily on volunteer anglers who were permitted to fish anywhere within Deadwater Slough. During a sampling occassion (day), anglers would boat or hike their catch (periodically or upon filling a livewell) to a processing station at the boat ramp approximately 500 m downstream from the top of the slough. For each fish, we recorded the date of capture, species, total length (TL; mm), and whether the fish was previously marked. Unmarked fish were given a physical mark (e.g., hole punch of lower caudal, upper caudal, left pelvic, right pelvic) unique to each day and released. Fish that died prior to release were still counted as captured (or recaptured), but not included in the number of marked fish available for recapture after that date. Released individuals were distributed throughout Deadwater Slough to facilitate mixing back into the population. Finally, we recorded the angling start and end time for each crew (person or combinations of persons) to facilitate calculations of CPUE.

Four different mark-recapture estimators were employed to explore variation in abundance estimates among them. Estimators fell into two broad categories: single census and multiple census. For the single census estimators, we treated the first week as the first (mark) sample, and the following (second) week as the second (recapture) sample, pooling data within each of those weeks. Alternatively, the multiple census estimators treat each day as a sample and use information about the total marked fish from all previous samples to infer the total abundance.

For the single census estimators, we first used the Lincoln-Peterson estimator:

\[ \hat{N} = \frac{(M)(n)}{(m)} \] where \(M\) is the total number of fish marked and returned to the population during the first sample (week), \(n\) is the total number of fish caught in the second sample (week), and \(m\) is the number of marked fish caught during the second week. Because the Lincoln-Peterson estimator can be biased with small samples, we also employed the Chapman-modified Lincoln-Peterson estimator:

\[ \hat{N} = \frac{(M + 1)(n + 1)}{(m + 1)} - 1 \] For a multiple census estimator, we first utilized the Schnabel method:

\[ \hat{N} = \frac{\sum\limits_{i = 1}^k n_i M_i}{\left(\sum\limits_{i = 1}^k m_i \right) + 1} \]

where \(M\), \(n\), and \(m\) are indexed by each sample (day), \(i\). With small numbers of recaptures, it may not be appropriate to calculate the standard error of the Schnabel estimator but the 95% confidence intervals can be calculated using a Poisson approximation (Krebs 1999). The Schnabel estimator is essentially a weighted average of a series of Lincoln-Peterson estimators (with a Champman modification). Lastly, we explored one additional multiple-census estimator, the Schumacher-Eschmeyer, which is based on minimizing the weighted sum of squares between the proportion of marked individuals in the sample and the unknown proportion of marked individuals in the population. The Schumacher-Eschmeyer estimator is as follows:

\[ \hat{N} = \frac{\sum\limits_{i = 1}^k n_i M^2_i}{\sum\limits_{i = 1}^k m_i M_i} \]

The mark-recapture estimators allowed us to estimate the abundance of Northern Pikeminnow during the fall 2019 and 2020 sample efforts (the estimate is of the initial population size). However, to estimate the abundance of Northern Pikeminnow during spring 2021, which was not a mark-recapture design, we used the ratio of total CPUE in the fall efforts to the total CPUE in the spring efforts, then multiplied that ratio by the average abundance from the fall efforts. This approach assumes equal capture probabilities between fall and spring efforts.

All mark-recapture abundance estimators assumed: (1) the population is closed (no immigration, emigration, births or deaths during the sampling period), (2) all fish have equal chance of being caught in the second (and subsequent) sample(s), (3) marking a fish does not affect its chances of recapture, (4) no loss of marks, and (5) marks are not missed or mistaken.

Northern Pikeminnow in Deadwater Slough were not only abundant, but anglers reported that many were large-bodied. For that reason, we calculated the proportional stock density (PSD) for Northern Pikeminnow in Deadwater Slough, which is a measure of species size structure. PSD is a ratio, typically expressed as a percentage, between the number of “quality-sized” (larger) individuals and “stock-sized” individuals:

\[ PSD_{i} = 100 * \frac{FQ_{i}}{FS_{i}} \]

where \(FQ_{i}\) is the number of fish \(\ge\) quality-length for species \(i\), and \(FS_{i}\) is the number of fish \(\ge\) stock-length for species \(i\). Quality and stock size definitions vary by species. We calculated the PSD for Northern Pikeminnow in Deadwater Slough using 380 mm TL for quality-size and 250 mm TL for stock-size (Winther et al. 2020).

Stomach Contents

Gastric lavage (Foster 1977) was used to examine the stomach contents of Northern Pikeminnow for the presence of juvenile Chinook Salmon and other fishes (e.g., juvenile steelhead, juvenile Sockeye Salmon, Redside Shiner Richardsonius balteatus, etc.). Immediately following lavage, stomach contents of individuals were preserved with 99% isopropyl alcohol in whirl-paks for later analysis in a controlled environment. For each sample, wet weight (grams) was recorded for the total combined stomach content, including all non-fish items (e.g., macroinvertebrates, organic matter), as well as that portion consisting of fish parts or matter. Fish and fish remnants were identified to the lowest taxonomic unit using diagnostic bones, or were categorized as unknown (Hansel et al. 1988; Frost 2000). A subset of Northern Pikeminnow captures (~ 5%) were euthanized for dissection after gastric lavage to validate the efficacy of the methodology.

Fish Consumption Potential

To estimate the total consumption potential (i.e., the number or amount of juvenile Chinook Salmon consumed during the peaks of DSR [fall] and NRR [spring] emigrations) of Northern Pikeminnow in Deadwater Slough, we used the Fish Bioenergetics v4.0 application developed by Deslauriers et al. (2017) and applied in the R statistical software (R Core Team 2021). The daily rate of consumption in grams for an individual Northern Pikeminnow was estimated based on the following inputs: predator and prey energy densities, predator start and end weights, and water temperatures. Predator energy density for Northern Pikeminnow was fixed at 6,703 Joules(J)/g (Deslauriers et al. 2017). Prey energy densities were fixed at 21,500 J/g for juvenile Chinook Salmon (Moss et al. 2016). The predator start weight was calculated using a weight-length formula (Parker et al. 1995) and the average length of Northern Pikeminnow caught in Deadwater Slough during our study. The average length of Northern Pikeminnow caught during the fall and spring efforts were 352.9 and 393.7 mm, respectively, which calculated to average starting weights of 504.4 and 670.5 g. We assumed no growth in individual Northern Pikeminnow, and so predator start and end weights are equal (i.e., we assumed the population biomass was stable). Mean daily water temperatures were summarized from 15-minute interval temperature readings available March 3, 2013 to June 14, 2021 from a gage station approximately 22 river kilometers downstream of Deadwater Slough near Shoup, Idaho (U.S. Geological Survey 2016).

The largest uncertainty among model parameters was the proportion of the Northern Pikeminnow diet consisting of fish versus non-fish prey items (e.g., invertebrates). Therefore, we conducted nine model runs with varying proportions of fish in the diet ranging from 10-90%, in 10% increments. Energy densities for invertebrates was fixed at 3,000 J/g. Because we were unable to differentiate juvenile Chinook Salmon from other fish prey consumed during model runs, we assume all fish prey have the same energy densities as juvenile Chinook Salmon i.e., juvenile Chinook Salmon were used as surrogates for all fish prey items. Separate models were also run for a 77-day fall period from September 15 - November 30 and a 92-day period in the spring from March 1 - May 31 to coincide with peak emigrations of DSR and NRR juveniles from the Lemhi River, the largest Chinook Salmon population in the Upper Salmon MPG. This resulted in 18 separate model runs (9 diet scenarios and 2 seasons). Each model run provided an estimate of the total amount (grams) of fish consumed by an individual Northern Pikeminnow during the fall or spring time periods. To estimate the total biomass of fish potentially consumed by all Northern Pikeminnow in Deadwater Slough, we multiplied the grams of fish consumed by an individual Northern Pikeminnow by the estimated predator population size during each period (fall and spring).

Potential Impacts to Chinook Salmon Populations

Multiplying the grams of fish consumed by an individual Northern Pikeminnow by the Northern Pikeminnow population size provides an estimate of total fish consumption potential for the population. However, the bioenergetics model runs do not differentiate juvenile Chinook Salmon from other fish prey items. Therefore, to understand potential impacts to local Chinook Salmon populations, we need to parse out the portion of total fish consumed that is juvenile Chinook Salmon. Unfortunately, estimates of the proportion of the Northern Pikeminnow fish prey that are juvenile Chinook Salmon are sparse, as well as information on the prevalence of juvenile Chinook Salmon relative to other fish prey items. Therefore, we instead consider hypothetical scenarios where Chinook Salmon are 1%, 10%, 30%, and 50% of fish prey consumed, by weight, which provide estimates of the total biomass (grams) of juvenile Chinook Salmon consumed by Northern Pikeminnow during the DSR and NRR emigration periods.

The estimate of total juvenile Chinook Salmon biomass consumed was converted to the total number of fish consumed each period by dividing by the average weight of DSR (10.3 g) and NRR (10.9 g) emigrants in the Upper Salmon MPG. The average weights of DSR and NRR emigrants were taken from fish caught at seven rotary screw traps located throughout the Upper Salmon MPG upstream of Deadwater Slough during the fall and spring time periods.

Results

Northern Pikeminnow Population Size

Using hook-and-line angling, we caught 664, 797, and 202 Northern Pikeminnow during the fall 2019, fall 2020, and spring 2021 efforts, respectively, for a total of 1,663 fish. The CPUE, measured as number of Northern Pikeminnow caught per angler hour, was 1.84 for fall 2019, 1.23 for fall 2020, and 0.81 for spring 2021 (Figure 2). During the mark-recapture efforts, there were eight recaptures in fall 2019 which included seven unique individuals. There were an additional six recaptures in fall 2020 consisting of five unique individuals. Summaries of parameters for both the single census and multiple census mark-recapture estimators are provided in Table 1 and Table 2, respectively.

Mark-recapture abundance estimates for Northern Pikeminnow in Deadwater Slough during fall 2019 ranged from 13,298 to 20,615 and for fall 2020 ranged from 24,381 to 42,475 (Table 3). In both cases, the multiple census estimators were larger than the single census estimators. Confidence intervals for the fall 2020 effort were more varied, likely due to fewer recaptures during that year.

All mark-recapture estimators yielded estimates with overlapping confidence intervals. However, because our sampling design most closely matched a multiple census estimator (i.e., we released marked fish each day and recaptures from all previous days were recorded for each day) and required a small sample size modification, we used the Schnabel estimates of abundance for further analyses and results. Accordingly, the mean of Northern Pikeminnow abundance for the two fall sampling events was 27,874 (95% CI: 14,244 - 59,388). By multiplying the ratio (0.53) between the spring CPUE (0.81) and mean of the fall CPUEs (1.53), we estimated there to be 14,799 (95% CI: 7,562 - 31,530) Northern Pikeminnow in Deadwater Slough during spring 2021.

The PSD for Northern Pikeminnow in Deadwater Slough across all three efforts was 50% i.e., 50% of captured fish were over the quality size classification defined by Winther et al. (2020). For comparison, Winther et al. (2020) reported a PSD for Northern Pikeminnow captured below Bonneville Dam of 41% and 18% in Bonneville Reservoir, suggesting that a larger fraction of Northern Pikeminnow in Deadwater Slough are of a quality size relative to populations reported elsewhere in the Columbia River. In fact, the state record Northern Pikeminnow, which measured 639 mm (25.2 in) TL, was caught in Deadwater Slough during the fall 2020 effort.

Stomach Contents

We examined the stomach contents of Northern Pikeminnow from Deadwater Slough using gastric lavage. To validate the efficacy of the method, we euthanized 75 Northern Pikeminnow after gastric lavage was completed and removed the remaining food items via dissection. We found that gastric lavage successfully removed food items in all dissected individuals, indicating the method was effective. Of the 1,558 Northern Pikeminnow sampled using gastric lavage, we found contents in some form in 350 (22.5%) stomachs and confirmed fish or fish parts in 44 (12.6%) of those. Northern Pikeminnow captured during the spring 2021 effort had a higher proportion of diets with fish contents (%) than fish captured in the fall 2019 (%) and fall 2020 (%) (Table 4). For those individuals with stomach contents, the average wet weight of total contents was 0.98 g (median = 0.25 g; range = 0.01-17.39 g). Overall, fish or fish remnants were 11.7% of the overall stomach contents examined.

Fish Consumption Potential

We estimated an individual, average-size (352.9 mm) Northern Pikeminnow to consume 45.23 g (range: 22.84 - 50.73 g) of fish during the DSR emigration to maintain their body size, depending on the proportion of diet consisting of fish versus invertebrates (Figure 5). Similarly, we estimated an average-size (393.7 mm) Northern Pikeminnow to consume 44.29 g (range: 22.39 - 49.74 g) of fish during the NRR emigration (Figure 5). Given a mean Northern Pikeminnow population size of 27,874 from fall estimates and 14,798.67 for spring 2021 and using the median value from diet scenarios (i.e., 50% of their diet consists of fish), we’d estimate that the Northern Pikeminnow population can consume 1,260,741 g (range: 636,642.2 - 1,414,048 g) fish during the DSR emigration and an additional 655,433.1 g (range: 331,342.2 - 736,085.8 g) during the NRR emigration in Deadwater Slough. In total, the Northern Pikeminnow population can consume 1,916,174 g of fish during the DSR and NRR emigrations to maintain their population size if we assume 50% of their diet consists of fish.

Potential Impacts to Chinook Salmon Populations

  1. Start with total g of fish consumed by Northern Pikemnnow population, summing DSR and NRR time periods.
  2. Parse out portion of fish total g consumed that is juvenile Chinook Salmon (1%, 10%, 30%, 50%).
  3. Divide by average weight of juvenile Chinook Salmon (10.6 g), the average of DSR (10.3 g) and NRR (10.9 g) emigrants.
  4. Estimate of total juvenile Chinook Salmon consumed.
  5. Multiply by SAR to estimate “adult equivalents.”

Discussion

The mean population size of Northern Pikeminnow in the Deadwater Slough during the fall season was estimated to be 27,874 and during the spring 2021 to be 14,798.67. Those estimates translate to linear densities of 14,898 and 7,909 which is 21 and 11 times greater, respectively, than densities reported by Beamesderfer and Rieman (1991) for the John Day Reservoir and approximately five and three times greater than projections for the lower Columbia River (Beamesderfer et al. 1996). Although estimated densities for Deadwater Slough are high given the relative size (volume) of the Salmon River at the slough relative to the Columbia River, high CPUE combined with low recapture rates (one in every 119 fish were recaptures) suggest this area may support exceptionally high predator densities. Further, Beamesderfer and Rieman (1991) acknowledge that their methods were unsuitable for sampling offshore, noting that water velocity, depth, irregular bottom contours, and barge traffic made sampling offshore ineffective. Conversely, at Deadwater Slough the maximum depth was approximately 10-15 m with a relatively homogenous bottom contour, making angling an effective method throughout the entire reach. Our estimates of Northern Pikeminnow areal density are also higher than those reported elsewhere in the Columbia River basin. We estimate the density of Northern Pikeminnow in Deadwater Slough, as fish per hectare, to be 2,334 for fall and 1,239 for spring. For comparison, the only other available estimates of areal density, that the authors are aware of, are 4.4 fish per hectare for John Day Reservoir (Beamesderfer and Rieman 1991) and 15 fish per hectare in Lake Washington (Bartoo 1977). Given the dearth of current data available in the Columbia River basin, it is unclear how Northern Pikeminnow abundance estimates in this study compare to elsewhere in the Columbia River Basin. Widener et al. (2021) has identified a need for more current information on piscine predator populations in the Columbia River.

Mark-Recapture Models

All mark-recapture abundance estimators suggest a large number of Northern Pikeminnow occupying Deadwater Slough. Several assumptions in our estimators may have influenced the magnitude of our results. First, all models assumed a closed population, however, it is possible that immigration and emigration occurred. Nevertheless, our estimators still provide unbiased estimates of abundance assuming the immigration and emigration rates are equal between marked and unmarked fish. Given the size of the sampling area, and the limited length of the sampling period (two weeks), that assumption is likely met. The assumption of a closed population can also affect results if mortality occurs for some marked fish released back to the population. During our field processing, we did not observe mortalities or injury post release; however, the potential for mortality amongst marked and released fish cannot be discounted. Next, our estimates assumed that the second (and subsequent) sample(s) are random samples of the population. Angling methods may have a size selection bias, thereby limiting our abundance estimates to a size threshold susceptible to angling, and possibly biasing them low. We also consider the validity of our assumption of equal catchability of individual fish between sampling events to have introduced potential bias in our results. Unequal catchability for individual fish is difficult to assess, especially if being caught once influences the likelihood of recapture. During our study, we recaptured two individual fish during multiple surveys indicating that marked fish were still susceptible to angling. We further recognize that marked and released fish may not have mixed thoroughly back into the population prior to recapture events. Although we attempted to release marked fish evenly throughout Deadwater Slough, we anecdotally observed anglers congregating in particular areas of the slough to socialize or exploit the “good” fishing locales. If marked fish were not recently released or mixed back into those areas that would result in marked fish being less susceptible to capture. Violations of the assumption of equal catchability between marked and unmarked fish may lead to an overestimate of abundance of Northern Pikeminnow in our study if marked fish have lower catchability.

The spring abundance estimate was smaller than fall abundance estimates. Our approach assumed equal capture probabilities between the fall and spring sample efforts. Yet, we suspect that the capture probability in the spring could be lower as higher spring flows reduce Northern Pikeminnow densities biasing our spring abundance estimate low. Therefore, we believe our spring abundance estimate to be conservative.

Stomach Contents

We performed gastric lavage on nearly all Northern Pikeminnow collected during this study. In addition to juvenile Chinook Salmon, species included Redside Shiner, Largescale Sucker Catostomus macrocheilus, sculpin Cottus spp., and Mountain Whitefish Prosopium williamsoni; although decomposition from digestion rendered many fish contents unidentifiable. While this approach provided useful information on prey species, we have concerns about using angling and gastric lavage to quantify diet composition in Northern Pikeminnow. First, previous studies have indicated that diets of angled fish differ from those collected by other methods (e.g., electrofishing) (Hodgson and Cochran 1988), and further, angled fish are more likely to have empty stomachs (Jurajda et al. 2016), perhaps because they are hungry and actively searching for food. Second, unlike other predators such as Smallmouth Bass and Walleye that are commonly sampled using gastric lavage (Kamler and Pope 2001), Northern Pikeminnow are cyprinids and therefore lack a true stomach. Previous research indicates that performing gastric lavage on cyprinid species can be ineffective at quantifying diet composition (Hartleb and Moring 1995), or alternately, requires a modification in methodology to flush food items out the vent of the fish (Wasowicz and Valdez 1994). Our intent, rather, was to simply document predation on juvenile Chinook Salmon in Deadwater Slough, and given the potential ineffectiveness of gastric lavage using angling and for Northern Pikeminnow, we consider our diet composition estimates to be conservative estimates i.e., we believe that a larger proportion of the Northern Pikeminnow diet is fish and juvenile Chinook Salmon than our gastric results suggest.

Bioenergetics

An average-size Northern Pikeminnow from Deadwater Slough is estimated to consume 45.23 g (range: 22.84 - 50.73 g) and 44.29 g (range: 22.39 - 49.74 g) of fish during fall and spring, respectively, to maintain its body weight. This equates to 89.52 g of fish consumed during the DSR and NRR peak emigrations, assuming the median value among diet scenarios (i.e., if 50% of their diet consists of fish). Consumption ranged from 45.23 for a diet comprised of 10% fish to 100.47 for a diet comprised of 90% fish.

Our approach assumed that all available prey consumed by Northern Pikeminnow have energy densities equal to juvenile Chinook Salmon. Other salmonid including juvenile steelhead and Sockeye Salmon are likely also available prey in Deadwater Slough, especially during the spring migration. Additionally, hatchery smolt releases (Chinook Salmon, steelhead, and Sockeye Salmon) are also prevalent in the mainstem Salmon River during the spring migration. Other studies have used a similar approach where the energy densities of different salmonids were generalized using an estimate from a single species (Petersen and Ward 1999). However, our study is focused on potential impacts to natural-origin DSR and NRR Chinook Salmon, in which case, let’s consider a thought exercise where DSR and NRR Chinook Salmon represent only 10% and 25% of total fish prey consumed by Northern Pikemnnow during these periods. Under these scenarios, the average Northern Pikeminnow would consume 8.95 g (range 4.52 – 10.05 g) if 10% of prey were natural-origin Chinook Salmon and 22.38 g (range 11.31 – 25.12) if 25%.

Our estimated consumption of Chinook Salmon juveniles by Northern Pikeminnow was calculated using a static predator population. If we assume an average juvenile Upper Salmon MPG Chinook Salmon of 10.6 g, this equates to 0.84 juvenile Chinook Salmon (range 0.43 – 0.94) under the 10% scenario and 2.11 juvenile Chinook Salmon (range 1.07 – 2.37) under the 25% scenario. Let’s now consider a stable Northern Pikeminnow population size of 10,000 (roughly the mean of fall and spring abundance estimates, divided by two to account for half of marked individuals not being available for recapture). In this case, the Northern Pikeminnow population in Deadwater Slough would be estimated to consume 8,400 (range 4,300 – 9,400) or 21,100 (range 10,700 – 23,700) under the 10% and 25% scenarios, respectively.

Impacts to Adult Returns

Our primary objective was to assess the impact of predation on juvenile Chinook Salmon in the Upper Salmon MPG. However, there is also interest in estimating the number of additional adult Chinook Salmon that might be expected to return to Lower Granite Dam if predation in the Deadwater Slough was reduced or eliminated. Our approach here was to multiply the total estimated juvenile Chinook Salmon consumed by an estimate of Granite-to-Granite smolt-to-adult survival rate (SAR) from McCann et al. (2019) for Upper Salmon River MPG Chinook Salmon. Assuming an SAR of 0.00614 (SE: 0.00051) we’ll consider the “adult equivalents” to our estimates of juveniles consumed. Under this scenario, the consumed juveniles would be equivalent to 51.6 adults (range 26.4 – 57.7) and 129.6 adults (range 65.7 – 145.6). However, we believe these 10% and 25% scenarios to be conservative, especially during the fall DSR emigration period. During spring, natural-origin steelhead and Sockeye Salmon are also actively emigrating from the Upper Salmon River and hatchery-origin releases of all three species (Chinook Salmon, steelhead, and Sockeye Salmon) are present in the river, and so 10% or 25% may be reasonable. However, during fall, DSR Chinook Salmon are likely the dominant prey species available to Chinook Salmon in Deadwater Slough as 1) no (or at least minimal) hatchery releases are present in the river and 2) fall emigrations of steelhead and Sockeye Salmon are less prominent than in Chinook Salmon. Note that this thought exercise only includes the Northern Pikeminnow population in Deadwater Slough even though they are present throughout the upper Salmon River, although likely at lower densities than in Deadwater Slough, and also ignores potential impacts to ESA-listed steelhead and Sockeye Salmon that are also a portion of the Northern Pikeminnow diet.

Avian Predation

Although not formally assessed in this study, avian predators including Great Blue Herons Ardea herodias and Bald Eagles Haliaeetus leucocephalus are another potential source of mortality for juvenile Chinook Salmon and other ESA-listed salmonids in the Deadwater Slough. The Deadwater Slough is recognized as an important bird watching and nesting area due to the associated upstream braided channels and riparian and backwater habitats (Deadwater Slough - Audubon Important Bird Areas). Several piscivorous bird species have been documented using Deadwater Slough including the Common Mergus merganser and Hooded Lophodytes cucullatus mergansers, the Great Blue Heron, the Double-crested Cormorant Phalacrocorax auritus, and the Belted Kingfisher Megaceryle alcyon (eBird 2021). During the initial sampling effort in fall 2019, a two-person crew walked the entire reach including surrounding and upstream areas scanning for passive integrated transponder (PIT) tags that are implanted into thousands of natural- and hatchery-origin salmonids in the Upper Salmon River, annually. During that informal survey, 10 PIT tags were recovered near active bird nests and in an upstream anastomizing reach where herons and eagles are prevalent, suggesting that mortality may have been a result of avian predation. The PIT tag histories in PTAGIS indicate these tags were implanted into a combination of juvenile Chinook Salmon (3), Sockeye Salmon (4), and steelhead (3). Two of the Chinook Salmon and one steelhead were natural-origin, the remaining were of hatchery-origin. Avian predation contributes a major component of the total mortality for yearling Chinook Salmon in some locations in the lower Snake and Columbia rivers, particularly at hydroelectric dams and within reservoirs (Evans et al. 2012; 2016); however, we did not observe large colonies of piscivorous birds within the study area. Although we did observe Belted Kingfishers, mergansers, and Great Blue Herons at Deadwater Slough during field efforts; Great Blue Herons were especially prevalent in the anastomizing reach immediately upriver of the slackwater reach. Additionally, there is documentation of individual Double-crested Cormorants (eBird 2021) at the Deadwater Slough; however, the site is not within their breeding range, rather, it is part of a migration corridor. Given the current avian species known to occupy Deadwater Slough, it is unlikely that avian predation on juvenile salmonids is comparable to elsewhere in the Columbia River basin with large piscivorous bird colonies. Nevertheless, we hypothesize that the reservoir-like conditions and associated habitat features at the Deadwater Slough may increase the probability of avian predation on juvenile Chinook Salmon and other ESA-listed salmonids from the many piscivorous birds known to use the site. Future studies of predation would benefit from estimates from piscivorous avian predators.

Management Implications

Recovery of the Upper Salmon MPG is desired to support local fisheries and economies. At least five of the eight populations must meet criteria set forth by McElhany et al. (2000) and Interior Columbia Technical Recovery Team (2007) for the MPG to be considered viable and for recovery of the Snake River Evolutionary Significant Unit. Spawner abundance is perhaps the most important metric considered in determining a population’s viability. We estimated that consumption of juvenile Chinook Salmon by Northern Pikeminnow in the Deadwater Slough potentially reduces natural-origin returns to upriver populations by 26.4 to 145.6 adults, annually. This is between 4.7% and 26.1% of the total mean adult escapement reported for the Upper Salmon River MPG in 2017 - 2019 (Kinzer et al. 2020). These estimates only take into account predation on juvenile Chinook Salmon by Northern Pikeminnow, but presumably, juvenile steelhead and Sockeye Salmon are also prey items in Deadwater Slough, including hatchery releases. Therefore, adult returns of steelhead and Sockeye Salmon are affected as well, including for hathchery populations which provide for recreational fishing opportunities. Consequently, reducing predation mortality at Deadwater Slough could potentially benefit multiple upriver natural and hatchery populations, including multiple ESA-listed species. Because the Deadwater Slough is part of the migratory pathway for emigrating salmonids, including for multiple species and populations, the impact of that benefit could be higher than tributary rehabilitation actions, which typically benefit a single population. Moreover, the deepened, slack water conditions that favor Northern Pikeminnow at Deadwater Slough are a result of the failure of a manmade mining reservoir dam. Taken together, the Deadwater Slough is a favorable candidate for management or restoration actions to benefit local Chinook Salmon populations.

We envision two potential management actions at reducing predation at Deadwater Slough: 1) removing the Dump Creek alluvial fan and 2) a local Northern Pikeminnow bounty program to encourage harvest in Deadwater Slough aimed at reducing the predator population size. Both actions have evident pros and cons. Removing the Dump Creek alluvial fan has the benefit of restoring the Salmon River to a more natural flow (e.g., increased water velocity, decreased depth) and bathymetry (e.g., heterogeneous substrate) that likely existed in the reach prior to the failure of the mining reservoir dam in Dump Creek. Restoring flow would presumably reduce Northern Pikeminnow densities to levels similar to upstream and downstream reaches where recent studies identified higher survival and transition probabilities than those at Deadwater Slough (Axel et al. (2015); Ackerman et al. (2018); Porter et al. (2019)). However, we do not know the feasibility of removing or alleviating the Dump Creek alluvial fan, and if it could be removed or altered, future disturbance events could potentially move additional loose sediment in the Dump Creek drainage, thereby creating a new alluvial fan. Additionally, restoring the flow would likely diminish the value of the reach as an important bird watching and nesting area.

Northern Pikeminnow bounty programs have proven successful at reducing population sizes and impacts on emigrating salmonids elsewhere (citation). A local bounty program in the Salmon River has the potential to boost the local economy (to an unknown degree) which has been negatively impacted by low returns of Chinook Salmon, steelhead, and Sockeye Salmon in recent decades. Many local businesses rely on license and gear sales for local sport fisheries, but low salmon and steelhead returns have eliminated or reduced the duration of fisheries in some years, and a local Northern Pikeminnow bounty program could allow the local economy to recoup some of those lost revenues. On the downside, bounty programs require continual effort and harvest to suppress the population size and so would require year of continued funding and local buy-in to remain effective. If funding or effort was alleviated, Northern Pikeminnow would likely recover to a higher population size without additional alterations to the habitat.

Conclusions

The failure of a mining reservoir dam in the late 19th century created an alluvial fan at the confluence of the Salmon River and Dump Creek which has backed up the Salmon River and created an almost 2-kilometer section of unnaturally slow and deep river, resulting in favorable conditions for Northern Pikeminnow. We estimate that greater than 10,000 Northern Pikeminnow occupy Deadwater Slough with densities greater than reported elsewhere. Additionally, the slow water velocity and lack of cover have created conditions that favor predation on juvenile salmonids emigrating through Deadwater Slough including Chinook Salmon i.e., juvenile Chinook Salmon have slower movement rates and little refugia to avoid predation. We estimate that predation from Northern Pikeminnow in Deadwater Slough may reduce adult returns by as much as 146 or 26% or recent returns to the Upper Salmon River MPG. Therefore, predation in Deadwater Slough likely has a consequential impact on ESA-listed Chinook Salmon populations in the Upper Salmon MPG.

Acknowledgements

The authors extend much appreciation to the many volunteers that assisted with field efforts including collaborators from the Bureau of Reclamation, Idaho Governor’s Office of Species Conservation, Idaho Department of Fish and Game, and the Lemhi Regional Land Trust, among others. Special thanks to Jared Barker, Brian Hamilton, Tulley Mackey, and Chelsea Welke for their help in the field. We further appreciate the administrative support and guidance from Caitlin Alcott and Inter-Fluve. This manuscript benefitted from a review by Sean Gibbs. Funding for this study was provided by the Bureau of Reclamation, Pacific Northwest Regional Office (contract No. 140R1021F0018) with assistance from the Idaho Governor’s Office of Species Conservation.

Literature Cited

Ackerman, M. W., G. A. Axel, R. A. Carmichael, and K. See. 2018. Movement and distribution of sp/sum Chinook Salmon presmolts in the mainstem Salmon River, pilot study. Technical report prepared for Idaho Governor’s Office of Species Conservation.
Axel, G. A., M. Peterson, C. C. Kozfkay, B. P. Sandford, M. G. Nesbit, B. J. Burke, K. E. Frick, and J. J. Lamb. 2015. Characterizing migration and survival between the upper Salmon River basin and Lower Granite Dam for juvenile Snake River sockeye salmon, 2014. Report of research by Fish Ecology Division, Northwest Fisheries Science Center, National Marine Fisheries Service, National Oceanic and Atmospheric Administration and Idaho Department of Fish and Game. Prepared for Division of Fish and Wildlife, Bonneville Power Administration. 36 pp.
Bartoo, N. W. 1977. Population parameter estimates and energy budgets for peamouth, northern squawfish, and yellow perch in Lake Washington. Doctoral Dissertation. University of Washington, Seattle.
Beamesderfer, C. P., D. L. Ward, and A. A. Nigro. 1996. Evaluation of the biological basis for a predator control program on Northern Pikeminnow (Ptychocheilus Oregonensis) in the Columbia and Snake rivers. Canadian Journal of Fisheries and Aquatic Sciences 53(12):2898–2908.
Beamesderfer, R. C., and B. E. Rieman. 1991. Abundance and distribution of Northern Squawfish, Walleyes, and Smallmouth Bass in John Day Reservoir, Columbia River. Transactions of the American Fisheries Society 120:439–447.
Copeland, T., D. A. Venditti, and B. R. Barnett. 2014. The importance of juvenile migration tactics to adult recruitment in stream-type Chinook Salmon populations. Transactions of the American Fisheries Society 143(6):1460–1475.
Deslauriers, D., S. R. Chipps, J. E. Breck, J. A. Rice, and C. P. Madenjian. 2017. Fish Bioenergetics 4.0: An R-based modeling application. Fisheries 42(11):586–596.
eBird. 2021. eBird: An online database of bird distribution and abundance [web application]. eBird, Cornell Lab of Ornithology, Ithaca, New York. Available: http://www.ebird.org. Accessed: November 10, 2021.
Emerson, J. L. 1973. A history of the Salmon National Forest. United States Forest Service. 194 pp.
Evans, A. F., N. J. Hostetter, D. D. Roby, K. Collis, D. E. Lyons, B. P. Sandford, R. D. Ledgerwood, and S. Sebring. 2012. Systemwide evaluation of avian predation on juvenile salmonids from the Columbia River based on recoveries of passive integrated transponder tags. Transactions of the American Fisheries Society 141(4):975–989.
Evans, A. F., Q. Payton, A. Turecek, B. Cramer, K. Collis, D. D. Roby, P. J. Loschl, L. Sullivan, J. Skalski, M. Weiland, and C. Dotson. 2016. Avian predation on juvenile salmonids: Spatial and temporal analysis based on acoustic and passive integrated transponder tags. Transactions of the American Fisheries Society 145(4):860–877.
Foster, J. R. 1977. Pulsed gastric lavage: An efficient method of removing the stomach contents of live fish. The Progressive Fish-Culturist 39(4):166–169.
Fresh, K. L., S. L. Schroder, and Mark I. Carr. 2003. Predation by Northern Pikeminnow on hatchery and wild Coho Salmon smolts in the Chehalis River, Washington. North American Journal of Fisheries Management 23(4):1257–1264.
Frost, C. N. 2000. A key for identifying preyfish in the columbia river based on diagnostic bones. U.S. Geological Survey, Western Fisheries Research Center, Columbia River Research Laboratory.
Hansel, H. C., S. D. Duke, P. T. Lofy, and G. A. Gray. 1988. Use of diagnostic bones to identify and estimate original lengths of ingested prey fishes. Transactions of the American Fisheries Society 117(1):55–62.
Hartleb, C. F., and J. R. Moring. 1995. An improved gastric lavage device for removing stomach contents from live fish. Fisheries Research 24(3):261–265.
Hodgson, J. R., and P. A. Cochran. 1988. The effect of sampling methodology on diet analysis in largemouth bass (Micropterus Salmoides). SIL Proceedings, 1922-2010 23(3):1670–1675.
Interior Columbia Technical Recovery Team. 2007. Viability criteria for application to interior Columbia basin salmonid ESUs. National Marine Fisheries Service, Northwest Fisheries Science Center, National Oceanographic and Atmospheric Administration.
Jurajda, P., K. Roche, I. Sedlacek, and L. Vsetickova. 2016. Assemblage characteristics and diet of fish in the shallow coastal waters of James Ross Island, Antarctica. Polar Biology 39(12):2299–2309.
Kamler, J. F., and K. L. Pope. 2001. Nonlethal methods of examining fish stomach contents. Reviews in Fisheries Science 9(1):1–11.
Kinzer, R., R. Orme, M. Campbell, J. Hargrove, and K. See. 2020. Report to NOAA Fisheries for 5-Year ESA Status Review: Snake River Basin Steelhead and Chinook Salmon population abundance, life History, and diversity metrics calculated from in-stream PIT-tag observations (SY2010-SY2019). IPTDSW (In-stream PIT-tag detection systems workgroup). 118 pp.
Knutsen, C. J., and D. L. Ward. 1999. Biological characteristics of Northern Pikeminnow in the lower Columbia and Snake rivers before and after sustained exploitation. Transactions of the American Fisheries Society 128(6):1008–1019.
Krebs, C. J. 1999. Ecological methodology. Addison Welsey Educational Publishers, Inc., Menlo Park, CA.
McCann, J., B. Chockley, E. Cooper, B. Hsu, G. Scheer, S. Haeseker, R. Lessard, T. Copeland, E. Tinus, A. Storch, D. Rawding, and M. DeHart. 2019. Comparative survival study of PIT-tagged spring/summer/fall Chinook, summer steelhead, and Sockeye. Draft 2019 annual report. BPA Project #19960200. Available at: https://www.fpc.org/documents/CSS/2019CSS_FullDRAFT.pdf.
McElhany, P., M. H. Ruckelshaus, M. J. Ford, T. C. Wainwright, and E. P. Bjorkstedt. 2000. Viable salmonid populations and the recovery of evolutionarily significant units. U.S. Department of Commerce, National Oceanic and Atmospheric Administration, Technical Memo. NMFS-NWFSC-42.:156 pp.
Moss, J., J. Murphy, E. Fergusson, and R. Heintz. 2016. Allometric relationships between body size and energy density of juvenile Chinook (Oncorhynchus Tshawytscha) and Chum (O. Keta) Salmon across a latitudinal gradient. North Pacific Anadromous Fish Commission Bulletin 6(1):161–168.
Murphy, C. A., J. D. Romer, K. Stertz, I. Arismendi, R. Emig, F. Monzyk, and S. L. Johnson. 2021. Damming salmon fry: Evidence for predation by non-native warmwater fishes in reservoirs. Ecosphere 12(9).
Nemeth, D. J., and R. B. Kiefer. 1999. Snake River spring and summer Chinook Salmon choice for recovery. Fisheries 24(10):16–23.
NOAA. 2017. ESA recovery plan for Snake River spring/summer Chinook Salmon (Oncorhynchus Tshawytscha) and Snake River basin steelhead (Oncorhynchus Mykiss).
Parker, R. M., M. P. Zimmerman, and D. L. Ward. 1995. Variability in biological characteristics of Northern Pikeminnow in the lower Columbia and Snake rivers. Transactions of the American Fisheries Society 124(3):335–346.
Petersen, J. H. 1994. Importance of spatial pattern in estimating predation on juvenile salmonids in the Columbia River. Transactions of the American Fisheries Society 123(6):924–930.
Petersen, J. H., and D. L. Ward. 1999. Development and corroboration of a bioenergetics model for Northern Pikeminnow feeding on juvenile salmonids in the Columbia River. Transactions of the American Fisheries Society 128(5):784–801.
Porter, N. J., M. W. Ackerman, T. Mackey, G. A. Axel, and K. E. See. 2019. Movement and distribution of Chinook Salmon presmolts in the mainstem Salmon River, 2018/2019 annual report. Prepared by Biomark, Inc. - Applied Biological Services and Fish Ecology Division, Northwest Fisheries Science Center, National Marine Fisheries Service, National Oceanic and Atmospheric Administration. 47 pp.
R Core Team. 2021. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria.
Roni, P., P. J. Anders, T. J. Beechie, and D. J. Kaplowe. 2018. Review of tools for identifying, planning, and implementing habitat restoration for Pacific salmon and steelhead. North American Journal of Fisheries Management 38(2):355–376.
Rubenson, E. S., D. J. Lawrence, and J. D. Olden. 2020. Threats to rearing juvenile Chinook Salmon from nonnative Smallmouth Bass inferred from stable isotope and fatty acid biomarkers. Transactions of the American Fisheries Society 149(3):350–363.
Shively, R. S., T. P. Poe, M. B. Sheer, and R. Peters. 1996. Criteria for reducing predation by Northern Squawfish near juvenile salmonid bypass outfalls at Columbia River dams. Regulated Rivers: Research & Management 12(4-5):493–500.
Tabor, R. A., R. S. Shively, and T. P. Poe. 1993. Predation of juvenile salmonids by Smallmouth Bass and Northern Pikeminnow in the Columbia River near Richland, Washington. North American Journal of Fisheries Management 13(4):831–838.
U.S. Geological Survey. 2016. National Water Information System data available on the world wide web (USGS Water Data for the Nation), acessed June 14, 2021, at URL https://waterdata.usgs.gov/nwis/uv?site_no=13307000.
Venditti, D. A., D. W. Rondorf, and J. M. Kraut. 2000. Migratory behavior and forebay delay of radio-tagged juvenile fall Chinook Salmon in a lower Snake River impoundment. North American Journal of Fisheries Management 20(1):41–52.
Ward, D. L., J. H. Petersen, and J. J. Loch. 1995. Index of predation on juvenile salmonids by Northern Squawfish in the lower and middle Columbia River and in the lower Snake River. Transactions of the American Fisheries Society 124:321–334.
Wasowicz, A., and R. A. Valdez. 1994. A nonlethal technique to recover gut contents of Roundtail Chub. North American Journal of Fisheries Management 14(3):656–658.
Widener, D. L., J. R. Faulkner, S. G. Smith, T. M. Marsh, and R. W. Zabel. 2021. Survival estimates for the passage of spring-migrating juvenile salmonids through Snake and Columbia River dams and reservoirs, 2020. U.S. Department of Commerce, NOAA Contract Report NMFS-NWFSC-CR-2021-05.
Winther, E., C. M. Barr, C. Miller, and C. Wheaton. 2020. Report on the predation index, predator control fisheries, and program evaluation for the Columbia River basin Northern Pikeminnow sport reward program, 2020 annual report. Pacific States Marine Fisheries Commission, Washington Department of Fish and Wildlife, Oregon Department of Fish and Wildife, Funded by: U.S. Deparment of Energy, Bonneville Power Administration, Environment, Fish, and Wildife. Project Number 1990-077-00. Contract Number 71866.
Wydoski, R. S., and R. R. Whitney. 2003. Inland fishes of Washington: Second edition, revised and expanded. University of Washington Press.
Zimmerman, M. P., and D. L. Ward. 1999. Index of predation on juvenile salmonids by Northern Pikeminnow in the lower Columbia River basin from 1994-96. Transactions of the American Fisheries Society 128(6):995–1007.

Tables

Table 1: Summary of Northern Pikeminnow angling survey data for single-census mark-recapture estimators. M is the total number of fish caught in the first sample (week), marked and returned to the population; n is the total number of fish caught in the second sample (week), and m is the number of previously marked fish caught in the second week.
Effort M n m
Fall 2019 267 396 7
Fall 2020 500 291 5
Table 2: Summary of Northern Pikeminnow angling survey data for multiple-census mark-recapture estimators. M is the total number of fish caught in previous samples (days), marked and returned to the population; n is the total number of fish caught in the current sample (day); m is the number of marked fish caught in that sample; u is the number of unmarked Northern Pikeminnow caught, and R is the number of marked fish returned to the population. Dates of sampling are also provided.
Effort Date M n m u R
Fall 2019 2019-11-12 0 29 0 29 28
Fall 2019 2019-11-13 28 146 0 146 146
Fall 2019 2019-11-14 174 93 1 92 93
Fall 2019 2019-11-19 266 149 2 147 132
Fall 2019 2019-11-20 396 104 1 103 77
Fall 2019 2019-11-21 472 143 4 139 118
Fall 2020 2020-10-20 0 173 0 173 170
Fall 2020 2020-10-21 170 188 1 187 187
Fall 2020 2020-10-22 356 104 0 104 102
Fall 2020 2020-10-23 458 41 0 41 41
Fall 2020 2020-10-27 499 42 0 42 41
Fall 2020 2020-10-28 540 47 1 46 46
Fall 2020 2020-10-29 585 157 4 153 156
Fall 2020 2020-10-30 737 45 0 45 45
Table 3: Point estimates, standard errors, and 95% confidence intervals for each of the mark-recapture estimators and for the fall sampling efforts. Standard errors are not available for the multiple-census estimators.
Effort Estimator N SE 95% Confidence Interval
Fall 2019 Petersen 15,105 5,658 7,331 - 37,569
Fall 2019 Chapman 13,298 4,322 6,898 - 27,893
Fall 2019 Schnabel 18,732 NA 10,057 - 37,851
Fall 2019 Schumacher-Eschmeyer 20,615 NA 14,393 - 36,313
Fall 2020 Petersen 29,100 12,902 12,470 - 89,622
Fall 2020 Chapman 24,381 9,066 11,547 - 55,761
Fall 2020 Schnabel 37,016 NA 18,430 - 80,924
Fall 2020 Schumacher-Eschmeyer 42,475 NA 22,501 - 378,128
Table 4: Summary of gastric lavage results including the total number of fish captured and gastric lavaged. For Northern Pikeminnow that were gastric lavaged, we show the number containing stomach contents, and of those, the number containing fish contents.
Effort Captured Lavaged Stomach Contents Fish Contents
Fall 2019 664 660 57 12
Fall 2020 797 793 188 25
Spring 2021 202 105 105 7
Total 1,663 1,558 350 44

Figures

Map of the Deadwater Slough study area. The high-resolution orthoimage portion surrounding the Deadwater Slough was generated from aerial images taken by an unmanned aircraft system (drone). The red polygon indicates the reach characterized by unnaturally slow water velocities and a deepened channel. The location of the Dump Creek delta is indicated.

Figure 1: Map of the Deadwater Slough study area. The high-resolution orthoimage portion surrounding the Deadwater Slough was generated from aerial images taken by an unmanned aircraft system (drone). The red polygon indicates the reach characterized by unnaturally slow water velocities and a deepened channel. The location of the Dump Creek delta is indicated.

The number of Northern Pikeminnow caught per angler effort for each sample effort.

Figure 2: The number of Northern Pikeminnow caught per angler effort for each sample effort.

Estimates of abundance of Northern Pikeminnow using different mark-recapture estimators and for the fall sampling events. Single census estimators are shown as circles; multiple census estimators are shown as triangles. Note the y-axis is truncated at 100 and the upper confidence interval for the Schumacker-Eschmeyer estimator in 2020 extends beyond that.

Figure 3: Estimates of abundance of Northern Pikeminnow using different mark-recapture estimators and for the fall sampling events. Single census estimators are shown as circles; multiple census estimators are shown as triangles. Note the y-axis is truncated at 100 and the upper confidence interval for the Schumacker-Eschmeyer estimator in 2020 extends beyond that.

Length frequency histogram of Northern Pikeminnow caught using hook-and-line angling during the study.

Figure 4: Length frequency histogram of Northern Pikeminnow caught using hook-and-line angling during the study.

The cumulative amount of fish (g) consumed by an individual Northern Pikeminnow during the peak DSR (fall) and NRR (spring) Chinook Salmon emigrations (top panel) and corresponding daily mean water temperatures (bottom panel). The black line is the median among diet scenarios (i.e., where 50% of the Northern Pikeminnow diet consists of fish) and the gray area shows the range among model scenarios (i.e., 10 - 90% of Northern Pikeminnow diet consists of fish). Daily mean water temperatures were summarized from six years of data available from the Shoup gage station.

Figure 5: The cumulative amount of fish (g) consumed by an individual Northern Pikeminnow during the peak DSR (fall) and NRR (spring) Chinook Salmon emigrations (top panel) and corresponding daily mean water temperatures (bottom panel). The black line is the median among diet scenarios (i.e., where 50% of the Northern Pikeminnow diet consists of fish) and the gray area shows the range among model scenarios (i.e., 10 - 90% of Northern Pikeminnow diet consists of fish). Daily mean water temperatures were summarized from six years of data available from the Shoup gage station.

The estimated number of DSR (fall) and NRR (spring) emigrants consumed by Northern Pikeminnow in Deadwater Slough (top panel) across a range of diet scenarios where 10-90% of the diet consists of juvenile Chinook salmon. The bottom panel provides the adult equivalents assuming a Lower Granite Dam smolt-to-adult rate (SAR) from McCann et al. 2019.

Figure 6: The estimated number of DSR (fall) and NRR (spring) emigrants consumed by Northern Pikeminnow in Deadwater Slough (top panel) across a range of diet scenarios where 10-90% of the diet consists of juvenile Chinook salmon. The bottom panel provides the adult equivalents assuming a Lower Granite Dam smolt-to-adult rate (SAR) from McCann et al. 2019.

Colophon

This report was generated on 2022-01-10 12:11:13 using the following computational environment and dependencies:

#> - Session info ---------------------------------------------------------------
#>  setting  value
#>  version  R version 4.1.1 (2021-08-10)
#>  os       Windows 10 x64 (build 19042)
#>  system   x86_64, mingw32
#>  ui       RTerm
#>  language (EN)
#>  collate  English_United States.1252
#>  ctype    English_United States.1252
#>  tz       America/Los_Angeles
#>  date     2022-01-10
#>  pandoc   2.14.0.3 @ C:/Program Files/RStudio/bin/pandoc/ (via rmarkdown)
#> 
#> - Packages -------------------------------------------------------------------
#>  package     * version date (UTC) lib source
#>  assertthat    0.2.1   2019-03-21 [1] CRAN (R 4.1.1)
#>  backports     1.3.0   2021-10-27 [1] CRAN (R 4.1.1)
#>  bit           4.0.4   2020-08-04 [1] CRAN (R 4.1.1)
#>  bit64         4.0.5   2020-08-30 [1] CRAN (R 4.1.1)
#>  bookdown      0.24    2021-09-02 [1] CRAN (R 4.1.1)
#>  broom         0.7.10  2021-10-31 [1] CRAN (R 4.1.2)
#>  bslib         0.3.1   2021-10-06 [1] CRAN (R 4.1.1)
#>  cachem        1.0.6   2021-08-19 [1] CRAN (R 4.1.1)
#>  callr         3.7.0   2021-04-20 [1] CRAN (R 4.1.1)
#>  cellranger    1.1.0   2016-07-27 [1] CRAN (R 4.1.1)
#>  cli           3.1.0   2021-10-27 [1] CRAN (R 4.1.1)
#>  colorspace    2.0-2   2021-06-24 [1] CRAN (R 4.1.1)
#>  crayon        1.4.2   2021-10-29 [1] CRAN (R 4.1.1)
#>  data.table    1.14.2  2021-09-27 [1] CRAN (R 4.1.1)
#>  DBI           1.1.1   2021-01-15 [1] CRAN (R 4.1.1)
#>  dbplyr        2.1.1   2021-04-06 [1] CRAN (R 4.1.1)
#>  desc          1.4.0   2021-09-28 [1] CRAN (R 4.1.1)
#>  devtools      2.4.3   2021-11-30 [1] CRAN (R 4.1.2)
#>  digest        0.6.28  2021-09-23 [1] CRAN (R 4.1.1)
#>  dplyr       * 1.0.7   2021-06-18 [1] CRAN (R 4.1.1)
#>  ellipsis      0.3.2   2021-04-29 [1] CRAN (R 4.1.1)
#>  english     * 1.2-6   2021-08-21 [1] CRAN (R 4.1.1)
#>  evaluate      0.14    2019-05-28 [1] CRAN (R 4.1.1)
#>  fansi         0.5.0   2021-05-25 [1] CRAN (R 4.1.1)
#>  farver        2.1.0   2021-02-28 [1] CRAN (R 4.1.1)
#>  fastmap       1.1.0   2021-01-25 [1] CRAN (R 4.1.1)
#>  forcats     * 0.5.1   2021-01-27 [1] CRAN (R 4.1.1)
#>  fs            1.5.0   2020-07-31 [1] CRAN (R 4.1.1)
#>  FSA         * 0.9.1   2021-07-17 [1] CRAN (R 4.1.1)
#>  generics      0.1.1   2021-10-25 [1] CRAN (R 4.1.1)
#>  ggplot2     * 3.3.5   2021-06-25 [1] CRAN (R 4.1.1)
#>  glue          1.4.2   2020-08-27 [1] CRAN (R 4.1.1)
#>  gridExtra   * 2.3     2017-09-09 [1] CRAN (R 4.1.1)
#>  gtable        0.3.0   2019-03-25 [1] CRAN (R 4.1.1)
#>  haven         2.4.3   2021-08-04 [1] CRAN (R 4.1.1)
#>  here        * 1.0.1   2020-12-13 [1] CRAN (R 4.1.1)
#>  highr         0.9     2021-04-16 [1] CRAN (R 4.1.1)
#>  hms           1.1.1   2021-09-26 [1] CRAN (R 4.1.1)
#>  htmltools     0.5.2   2021-08-25 [1] CRAN (R 4.1.1)
#>  httr          1.4.2   2020-07-20 [1] CRAN (R 4.1.1)
#>  janitor     * 2.1.0   2021-01-05 [1] CRAN (R 4.1.1)
#>  jquerylib     0.1.4   2021-04-26 [1] CRAN (R 4.1.1)
#>  jsonlite      1.7.2   2020-12-09 [1] CRAN (R 4.1.1)
#>  kableExtra  * 1.3.4   2021-02-20 [1] CRAN (R 4.1.1)
#>  knitr         1.36    2021-09-29 [1] CRAN (R 4.1.1)
#>  labeling      0.4.2   2020-10-20 [1] CRAN (R 4.1.1)
#>  lattice       0.20-44 2021-05-02 [2] CRAN (R 4.1.1)
#>  lemon       * 0.4.5   2020-06-08 [1] CRAN (R 4.1.2)
#>  lifecycle     1.0.1   2021-09-24 [1] CRAN (R 4.1.1)
#>  lubridate   * 1.8.0   2021-10-07 [1] CRAN (R 4.1.1)
#>  magrittr    * 2.0.1   2020-11-17 [1] CRAN (R 4.1.1)
#>  memoise       2.0.1   2021-11-26 [1] CRAN (R 4.1.2)
#>  modelr        0.1.8   2020-05-19 [1] CRAN (R 4.1.1)
#>  munsell       0.5.0   2018-06-12 [1] CRAN (R 4.1.1)
#>  pillar        1.6.4   2021-10-18 [1] CRAN (R 4.1.1)
#>  pkgbuild      1.2.1   2021-11-30 [1] CRAN (R 4.1.2)
#>  pkgconfig     2.0.3   2019-09-22 [1] CRAN (R 4.1.1)
#>  pkgload       1.2.4   2021-11-30 [1] CRAN (R 4.1.2)
#>  plyr          1.8.6   2020-03-03 [1] CRAN (R 4.1.1)
#>  prettyunits   1.1.1   2020-01-24 [1] CRAN (R 4.1.1)
#>  processx      3.5.2   2021-04-30 [1] CRAN (R 4.1.1)
#>  ps            1.6.0   2021-02-28 [1] CRAN (R 4.1.1)
#>  purrr       * 0.3.4   2020-04-17 [1] CRAN (R 4.1.1)
#>  R6            2.5.1   2021-08-19 [1] CRAN (R 4.1.1)
#>  Rcpp          1.0.7   2021-07-07 [1] CRAN (R 4.1.1)
#>  readr       * 2.0.2   2021-09-27 [1] CRAN (R 4.1.1)
#>  readxl        1.3.1   2019-03-13 [1] CRAN (R 4.1.1)
#>  remotes       2.4.2   2021-11-30 [1] CRAN (R 4.1.2)
#>  reprex        2.0.1   2021-08-05 [1] CRAN (R 4.1.1)
#>  rlang         0.4.12  2021-10-18 [1] CRAN (R 4.1.1)
#>  rmarkdown     2.11    2021-09-14 [1] CRAN (R 4.1.1)
#>  rprojroot     2.0.2   2020-11-15 [1] CRAN (R 4.1.1)
#>  rstudioapi    0.13    2020-11-12 [1] CRAN (R 4.1.1)
#>  rvest         1.0.2   2021-10-16 [1] CRAN (R 4.1.1)
#>  sass          0.4.0   2021-05-12 [1] CRAN (R 4.1.1)
#>  scales        1.1.1   2020-05-11 [1] CRAN (R 4.1.1)
#>  sessioninfo   1.2.2   2021-12-06 [1] CRAN (R 4.1.2)
#>  snakecase     0.11.0  2019-05-25 [1] CRAN (R 4.1.1)
#>  stringi       1.7.5   2021-10-04 [1] CRAN (R 4.1.1)
#>  stringr     * 1.4.0   2019-02-10 [1] CRAN (R 4.1.1)
#>  svglite       2.0.0   2021-02-20 [1] CRAN (R 4.1.1)
#>  systemfonts   1.0.3   2021-10-13 [1] CRAN (R 4.1.1)
#>  testthat      3.1.0   2021-10-04 [1] CRAN (R 4.1.1)
#>  tibble      * 3.1.5   2021-09-30 [1] CRAN (R 4.1.1)
#>  tidyr       * 1.1.4   2021-09-27 [1] CRAN (R 4.1.1)
#>  tidyselect    1.1.1   2021-04-30 [1] CRAN (R 4.1.1)
#>  tidytable   * 0.6.5   2021-09-03 [1] CRAN (R 4.1.1)
#>  tidyverse   * 1.3.1   2021-04-15 [1] CRAN (R 4.1.1)
#>  tzdb          0.2.0   2021-10-27 [1] CRAN (R 4.1.1)
#>  usethis       2.1.3   2021-10-27 [1] CRAN (R 4.1.1)
#>  utf8          1.2.2   2021-07-24 [1] CRAN (R 4.1.1)
#>  vctrs         0.3.8   2021-04-29 [1] CRAN (R 4.1.1)
#>  viridisLite   0.4.0   2021-04-13 [1] CRAN (R 4.1.1)
#>  vroom         1.5.5   2021-09-14 [1] CRAN (R 4.1.1)
#>  webshot       0.5.2   2019-11-22 [1] CRAN (R 4.1.1)
#>  withr         2.4.3   2021-11-30 [1] CRAN (R 4.1.2)
#>  xfun          0.27    2021-10-18 [1] CRAN (R 4.1.1)
#>  xml2          1.3.2   2020-04-23 [1] CRAN (R 4.1.1)
#>  yaml          2.2.1   2020-02-01 [1] CRAN (R 4.1.1)
#> 
#>  [1] C:/Users/seek1477/OneDrive - Washington State Executive Branch Agencies/Documents/R/win-library/4.1
#>  [2] C:/Program Files/R/R-4.1.1/library
#> 
#> ------------------------------------------------------------------------------

The current Git commit details are:

#> Local:    main C:/Users/seek1477/OneDrive - Washington State Executive Branch Agencies/Documents/Git/MyProjects/DeadwaterPaper
#> Remote:   main @ origin (https://github.com/mackerman44/DeadwaterPaper.git)
#> Head:     [4e80073] 2022-01-10: updated script to examine the impacts if not all marked fish were available for recapture.